Long standing problem of 210 Bi and the realistic neutron-proton effective interaction 
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The odd-odd nucleus 210 Bi is studied within the framework of the shell model using effective 
two-body matrix elements derived from the CD-Bonn nucleon-nucleon potential. The experimental 
energies of the proton-neutron multiplet 7r/ig/ 2 ^<?9/2 are remarkably well reproduced by the theory, 
which accounts for the 1~ state being the ground state instead of the 0~ predicted by the Nordheim 
strong coupling rule, ft is shown that the core-polarization effects are crucial to produce this 
inversion. The similarity between neutron-proton multiplets in the 132 Sn and 208 Pb regions is 
discussed in connection with the effective interaction. 

PACS numbers: 21.60.Cs, 21.30.Fe, 27.60.+W 



A fundamental goal of nuclear physics is to under- 
stand the properties of nuclei starting from the forces be- 
tween nucleons. Within the framework of the shell model, 
which is the basic framework for the study of complex 
nuclei in terms of nucleons, this implies the derivation 
of the two-body effective interaction V c g from the free 
nucleon-nucleon (NN) potential. This has long been a 
main problem of microscopic nuclear structure theory. 
Suffice to say that since the early 1950s through the mid 
1990s in the vast majority of shell-model calculations ei- 
ther empirical effective interactions containing adjustable 
parameters have been used or the two-body matrix ele- 
ments themselves have been treated as free parameters. 

Great progress toward "realistic" shell- model calcula- 
tions has been achieved over the last decade, proving the 
ability of effective interactions derived from modern NN 
potentials to provide, with no adjustable parameters, an 
accurate description of nuclear properties. This makes 
it very interesting to put to the test realistic shell-model 
calculations on peculiar nuclear properties which could 
not be unambiguously explained by earlier calculations. 
One such case is provided by the nucleus 210 Bi. 

The nucleus 210 Bi with one proton and one neutron 
outside doubly magic 208 Pb core has long been the sub- 
ject of theoretical studies. In this nucleus there are ten 
states below 600 keV with angular momenta from zero 
to nine, which in terms of the shell model constitute 
a neutron-proton multiplet arising from the 7r/ig/2i / <?9/2 
configuration. The experimental determination of the 
ground-state spin of 210 Bi, J 71 " = 1 _ , dates back to some 
fifty years ago [1| and since then it has been a challenge to 
shell-model theorists to explain this feature, which is at 
variance with the prediction, J 71 " = _ , of the Nordheim 
strong coupling rule [3. 

Early calculations [!, 0, [j| made use of empirical 
two-body interactions containing adjustable parameters. 
These studies led to the conclusion that tensor-force com- 
ponents were needed to explain the inversion of the 0~ 



and 1" states. Realistic shell-model calculations start- 
ing from the Hamada- Johnston NN potential were 
performed for 210 Bi by Herling and Kuo Q in the early 
1970s. The hard core of the Hamada-Johnston poten- 
tial was handled by using the Brueckner reaction matrix 
G. Then, the effective two-body interaction was con- 
structed by adding to the G matrix certain classes of 
diagrams arising from core-excitation processes, account- 
ing in this way for configurations left out of the space in 
which the shell-model calculation was performed. A good 
overall agreement between the calculated and experimen- 
tal spectra was obtained in [7], but the above inversion 
was not explained. Some twenty years later, calculations 
in the 208 Pb region were performed by Warburton and 
Brown with some modifications to the Kuo-Herling 
interaction. In particular, for the 7r/ig/ 2 ^39/2 multiplet, 
modifications of all the corresponding two-body matrix 
elements were introduced to reproduce the experimental 
energies. 

In the work of [9] the results of a shell-model calcu- 
lation for nuclei in the lead region including 210 Bi were 
presented together with a brief but comprehensive review 
of studies on this nucleus through 1972. The authors of 
Q made use of an empirical effective interaction with no 
tensor components but containing multipole forces with 
adjustable strengths to account for model-space trunca- 
tion effects. Their conclusion was that an unambiguous 
explanation of the inversion of the 0~ and 1" states still 
remained to be given. In fact, they showed that their 
model was not able to simultaneously describe the correct 
ordering of the two states and the energies of the higher- 
lying states of the ithg^vgajo multiplet. They pointed 
out that, while the studies evidenced the key role 

of the tensor force to reproduce the whole multiplet, real- 
istic effective interactions as that used in Q had failed to 
reproduce the 210 Bi ground state although the starting 
NN potential contained explicitly tensor components. 

This situation had not changed much after more than 
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two decades, as is reflected in the work of Ref. [l(| , where 
a description of 210 Bi is given by using a phenomenolog- 
ical model, which, besides an empirical neutron-proton 
effective interaction, includes explicitly macroscopic de- 
grees of freedom through a 208 Pb vibrating core. 

As mentioned above, in the last ten years or so shell- 
model calculations employing realistic effective interac- 
tions derived from modern NN potentials have produced 
results in remarkably good agreement with experiment 
for a number of nuclei in various mass regions. In the 
most recent calculations, the difficulty of dealing with po- 
tentials having a strong repulsive core has been overcome 
by resorting to the so called Mow-k approach Within 
this approach a low-momentum NN potential, Vjow-k, is 
obtained by integrating out the high-momentum modes 
of the NN potential down to a cutoff momentum A, with 
the requirement that the deuteron binding energy as well 
as the half-on-shell T matrix of the original potential are 
preserved. This T-matrix equivalence procedure, as de- 
scribed in detail in Ref. [ll|, produces a smooth poten- 
tial which is suitable for being used directly in nuclear 
structure calculations. The shell- model effective interac- 
tion is then constructed within the framework of a Q-box 
folded-diagrams expansion [l2| , namely it is expressed as 
a Q-box, which is vertex function composed of irreducible 
valence-linked diagrams at any order in Vi ow -k, plus the 
folded-diagram series. Note that the G-matrix calcula- 
tion of the effective interaction performed by Herling and 
Kuo [3] is equivalent to retain the first three terms of 
the Q-box without folded diagrams, i.e., the G matrix 
and the core-polarization diagrams corresponding to one 
particle-one hole (lplh) and two particle-two hole (2p2h) 
excitations. 

In recent years, particular attention from the exper- 
imental and theoretical point of view has been focused 
on nuclei around doubly magic 132 Sn. In this region, the 
counterpart of 210 Bi is the nucleus 134 Sb with one proton 
and one neutron in the Z = 50 — 82 and N = 82 — 126 
shells, respectively. In this case, the ground state has 
J n = 0~ and is nearly degenerate with the first-excited 
J n = 1~ state, the latter lying at 13 keV. These two 
states are members of the lowest neutron-proton mul- 
tiplet arising from the Trg^ / 2 v 'fi /2 configuration. Very 
recently, we have performed a shell-model study of this 
nucleus [l3[ making use of an effective interaction derived 
from the CD-Bonn NN potential [3] . Our results turned 
out to be in very good agreement with the experimental 
data including the very small energy spacing between the 
0~ and 1~ states. A main finding of this study was that 
core polarization effects, in particular those arising from 
lplh excitations, introduce modifications in the effective 
interaction which are essential to reduce the spacing be- 
tween the 1~ and CP states. 

The above achievement is at the origin of the present 
paper, where we present the results of a realistic shell- 
model calculation for 210 Bi conducted along the same 



lines as that for 134 Sb. Note that each of the two 
single-particle levels composing the ^57/2^/7/2 multiplet 
in 134 Sb has a counterpart with the same radial quantum 
number and one more unit in the angular momenta I and 
j in the nh 9 / 2 vgg/2 multiplet of 210 Bi. As we shall see, 
this multiplet is very well reproduced with the 0" and 
l - states in correct order. 

Our study of 210 Bi has been performed assuming 208 Pb 
as a closed core and taking as model space for the valence 
proton and neutron the six levels of the 82-126 shell and 
the seven levels of the 126-184 shell, respectively. 

The two-body effective interaction V e g has been de- 
rived by means of the Q-box folded-diagrams expansion 
(l2l | from the CD-Bonn NN potential, the short-range 
repulsion of the latter being renormalized by use of the 
Viow-k potential As in our previous studies for nu- 
clei around 132 Sn [H [H EH, the cutoff momentum A 
is given the value 2.2 fm _1 . This Viow-k) with addition 
of the Coulomb interaction for protons, is then used to 
calculate the Q-box . In our calculation, we include one- 
and two-body diagrams up to second order in the interac- 
tion, as they are explicitly shown in [l7j |. The computa- 
tion of these diagrams is performed within the harmonic- 
oscillator basis using intermediate states composed of all 
possible hole states and particle states restricted to the 
five proton and neutron shells above the Fermi surface. 
This guarantees the stability of the results when increas- 
ing the number of intermediate particle states. The os- 
cillator parameter is 6.88 MeV, as obtained from the ex- 
pression hu = 45 A^ 1 / 3 — 25 A -2 / 3 . Once the Q-box is 
calculated, the series of folded diagrams is summed up to 
all orders using the Lee-Suzuki iteration method [l8(. 

The effective interaction obtained by this procedure 
contains one- and two-body contributions, the former 
determining, once summed to the unperturbed Hamil- 
tonian, the single-particle energies. We have used the 
subtraction procedure of Ref. [17| so as to retain only 
the two-body term V e g, while the single-particle energies 
have been taken from experiment. The single- proton and 
-neutron energies adopted in our calculation are reported 
in Table I. They are relative to 208 Pb as obtained from 
the spectra of 209 Bi and 208 Pb [3| using mass excesses 
from fl]]. 

Let us now come to the results of our calculations, 
which have been performed by using the OXBASH shell- 
model code [2l[ . As for the binding energy of the ground 
state, our calculated value is 8.347 ± 0.003 MeV, which 
compares very well with the experimental one, 8.404 ± 
0.002 MeV 0. Note that the error on the calculated 
value arises from the experimental errors on the proton 
and neutron separation energies of 209 Bi and 209 Pb (io| . 

The calculated wave functions of the yrast states with 
J w from 0~ to 9~ are all dominated by the irhg/2 l/ 99/2 
configuration with very little configuration mixing. In 
fact, the percentage of this component ranges from 96% 
to 100% except for the 1 _ state where it has the minimum 
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TABLE I: Proton and neutron single-particle energies (in 
MeV). 



7r(n, l,j) 


6 


v(n, I, j) 


6 


O/lo/o 


-3.80 


lgg/2 


-3.94 


I/7/2 


-2.90 


0in/2 


-3.16 


0«13/2 


-2.19 


0jl5/2 


-2.51 


1/5/2 


-0.97 


2d 5 /2 


-2.37 


2P3/2 


-0.68 


3Sl/2 


-1.90 


2pi/ 2 


-0.16 


137/2 


-1.44 






2d 3 /2 


-1.40 



value of 91%. These states represent therefore the ten 
members of the 77/19/2^(79/2 multiplet and their excitation 
energies are reported in Fig. [TJ where they are compared 
with the experimental yrast states with the same spin 
and parity [l9(. Note that besides the ten members of 
the multiplet, only a second 1~ state is found below 600 
keV in both the experimental and calculated spectra. 
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FIG. 1: Proton-neutron 77/19/2^(79/2 multiplet in 210 Bi. The 
theoretical results are represented by open circles while the 
experimental data by solid triangles. 

From Fig. [TJ we see that the experimental pattern is 
well reproduced by the theory. We correctly predict the 
energy decrease from the 0~ to the 1~ state as well as the 
sharp and slight increase occurring in the next two steps. 
Then, from the 3 _ state on the experimental and the- 
oretical patterns stagger with the same magnitude and 
phase. The results shown in Fig. [T] are also in remark- 
ably quantitative agreement with the experimental data. 
In fact, the discrepancies are less than 50 keV for all the 
states except the 9~, which is predicted to lie at 90 keV 
above the experimental one. 

Let us now discuss our effective interaction. In our 
study of 134 Sb [l3| the effective two-body interaction 
was derived using precisely the same approach as that 
adopted here for 210 Bi and it turned out that the 
7r 37/2 l/ /7/2 multiplet was very well described. In this 
context, we analyzed the various terms contributing to 



the effective interaction in order to understand their rel- 
ative importance in determining the neutron-proton ma- 
trix elements leading to a successful description of the 
multiplet. This provided clear evidence of the crucial 
role played by the lplh excitations especially in regard 
to the 0~ — \~ spacing. 
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FIG. 2: Diagonal matrix elements of the two-body effective 
interaction for the 77/2,9/2^(79/2 configuration. 
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FIG. 3: Diagonal matrix elements of Vj ow _k and contributions 
from the two-body second-order diagrams for the 77/2.9/2^(79/2 
configuration. See text for comments. 

Here, we have performed the same analysis for the 
7r/ig/2JA7g/2 multiplet in 210 Bi. In Fig. [21 we report the 
diagonal matrix elements of the effective interaction as 
a function of J for the 77/19/2^(79/2 configuration. As ex- 
pected, we see that their behavior is quite similar to that 
shown in Fig. [JJfor the energies of the multiplet. Then, 
we show in Fig. [3] the four two-body second-order con- 
tributions to the effective interaction, namely the corre- 
sponding matrix elements of Vjow-k and the renormal- 
izations due to core polarization through lplh and lplh 
excitations as well as to the ladder diagrams which ac- 
count for excluded configurations above the chosen model 
space. As in the case of 134 Sb, the folding contribution 
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turns out to be quite irrelevant for the behavior of the 
neutron-proton multiplet. From Fig. [3]it is clear that the 
features of the neutron-proton effective interaction evi- 
denced in our study of 134 Sb persist in the 208 Pb region. 
In fact, we see that also in this region the neutron-proton 
matrix elements of Viow-k do not follow the experimental 
behavior, the most attractive one being that for = CP. 
This behavior is significantly modified by the Vi p ih term, 
which is always repulsive except for J 71 " = 1 — . 

It is worth mentioning that in a preliminary calcula- 
tion on 210 Bi [22J the quality of agreement between theory 
and experiment was not as good as that obtained here. 
In particular, we did not predict the CP and 1 _ levels in 
correct order. However, also in that case the contribution 
from lplh excitations was essential to reduce the spacing 
between the 1~ and CP states, although not sufficient to 
produce the inversion of these two states. It turned out 
that the correct ground state could be obtained increas- 
ing by a factor of about 2.5 the diagonal matrix elements 
of Vipih for the 7r/ig/2 ^179/2 configuration. 

In the present calculation, the lplh contribution brings 
about the right effect without the introduction of any 
modification. This is traced to the improvement in the 
calculation of the effective interaction resulting from the 
increased number of intermediate states used in the com- 
putation of the Q-box diagrams. The same result was 
also found for 134 Sb [H . 

At this point, we may conclude that our effective in- 
teraction contains the noncentral components needed to 
give the CP and l - levels in correct order. These arise 
from virtual interactions of the core particles induced by 
the NN potential and their weight with respect to the 
central components may be quite different from that rel- 
ative to the NN potential. 

In summary, we have shown that a shell-model inter- 
action derived from a modern NN potential without any 
modification leads to a direct solution of the long stand- 
ing problem of 210 Bi. 

This work was supported in part by the Italian 
Ministcro dell'Istruzione, dell'Universita e della Ricerca 
(MIUR). 



[1] J. R. Erskine, W. W. Buechner, and H. A. Enge, Phys 

Rev. 128, 720 (1962). 
[2] L. W. Nordheim, Phys Rev. 78, 294 (1950). 
[3] Y. E. Kim and J. 6. Rasmussen, Nucl. Phys. 47, 184 

(1963). 

[4] P. A. Mello and J. Flores, Nucl. Phys. 47, 177 (1963). 
[5] T. A. Hughes, R. Snow, and W. T. Pinkston, Nucl. Phys. 

82, 129 (1966). 
[6] T. Hamada and I. D. Johnston, Nucl. Phys. 34, 382 

(1962). 

[7] G. H. Herling and T. T. S. Kuo, Nucl. Phys. A 181, 113 
(1972). 

[8] E. K. Warburton and B. A. Brown, Phys. Rev. C 43, 602 
(1991). 

[9] C. W. Ma and W. W. True, Phys. Rev. C 8, 2313 (1973). 
[10] P. Alexa, J. Kvasil, N. V. Minh, and R. K. Sheline, Phys. 

Rev. C 55, 179 (1997). 
[11] S. Bogner, T. T. S. Kuo, L. Coraggio, A. Covello, and N. 

Itaco, Phys. Rev. C 65, 051301 (2002). 
[12] T. T. S. Kuo and E. Osnes, Lecture Notes in Physics, 

Vol. 364, (Springer- Verlag, Berlin, 1990). 
[13] L. Coraggio, A. Covello, A. Gargano, and N. Itaco, Phys. 

Rev. C 73, 031302(R) (2006). 
[14] R. Machleidt, Phys. Rev. C 63, 024001 (2001). 
[15] L. Coraggio, A. Covello, A. Gargano, and N. Itaco, Phys. 

Rev. C 72, 057302 (2005). 
[16] A. Covello, L. Coraggio, A. Gargano, and N. Itaco, Prog. 

Part. Nucl. Phys. 59, 401 (2007). 
[17] J. Shurpin, T. T. S. Kuo, and D. Strottman, Nucl. Phys. 

A 408, 310 (1983). 
[18] K. Suzuki and S. Y. Lee, Prog. Theor. Phys. 64,2091 

(1980). 

[19] Data extracted using the NNDC On-line Data Service 

from the XUNDL database, file revised as of July 6, 2007. 
[20] G. Audi, A. H. Wapstra, and C. Thibault, Nucl. Phys. A 

729, 337 (2003). 
[21] B. A. Brown, A. Etchegoyen, and W. D. M. Rae, The 

computer code OXBASH, MSU-NSCL, Report No. 524. 
[22] A. Covello, L. Coraggio, A. Gargano, and N. Itaco, Phys. 

Atom. Nucl. 67, 1611 (2004). 



